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Abstract. The Alpha Magnetic Spectrometer (AMS) 
is a particle physics detector designed for a high 
precision measurement of cosmic rays in space. AMS 
phase 2 (AMS-02) is scheduled to be installed on the 
ISS for at least three years from September 2010. The 
AMS-01 precursor experiment operated successfully 
during a 10-day NASA shuttle flight in June 1998. 
The orbital inclination was 51.7° at a geodetic al- 
titude between 320 to 380 km. Nearly 200,000 Z>2 
nuclei were observed by AMS-01 in the rigidity range 
1-40 GV. Using these data, it is possible to investigate 
the relative abundances and the energy spectra of 
the primary cosmic rays, providing relations with 
their sources and propagation processes. Preliminary 
results on the B/C ratio in 0.4-19 GeV/nucleon kinetic 
energy will be presented. 

Keywords: AMS-01, Spectrometer, Nuclei 



I. Introduction 

The boron to carbon nuclear ratio is very important for 
understanding the cosmic ray propagation in the Galaxy. 
Boron and carbon are respectively a pure secondary 
and a pure primary cosmic ray species. Since B is 
directly produced by C spallation in the interstellar 
medium, their relative abundance is sensitive to the 
traversed matter from the source, and the B/C energy 
distribution is an interesting quantity to study cosmic 
rays diffusion and interactions in the galactic medium. 
A precise measurement of the B/C ratio constraints the 
propagation parameters and, as a consequence, provides 
a better understanding of the astrophysical processes 
at the sources. The AMS-01 measurement of the B/C 
ratio is presented here in the kinetic energy range 0.4- 
19 GeV/n. 

II. The AMS-01 detector 

The AMS-01 spectrometer was composed of a cylin- 
drical permanent magnet, a silicon microstrip tracker, 
time of flight (TOF) scintillator planes, a Cerenkov 
counter and anti-coincidence counters [lj. The magnet 
(inner diameter 1.1m) provided a central dipole field 
with an analyzing power BL 2 — 0.14Tm 2 . Six layers 
of double-sided silicon microstrip tracker measured the 
trajectory of charged particles with an accuracy of 10/j,to 
(30/im) in the (non)bending coordinates, as well as 
providing multiple energy loss measurements. The TOF 
system had two orthogonal planes at each end of the 



magnet, and measured the particle transit time with an 
accuracy of ^120 psec; they also yielded energy loss 
measurement up to \Z\=2. A layer of anti-coincidence 
scintillation counters lined the inner surface of the 
magnet. A thin carbon layer were used as a shield to 
absorb low energy particles. 
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Fig. 1: The AMS-01 detector 

The detector response was simulated under a 
GEANT3-based software |2j. The effects of energy loss, 
multiple scattering, nuclear interactions and decays were 
included in the code, as well as efficiencies, resolutions 
and reconstruction algorithms. 

III. Event reconstruction and selection 

The main physical characteristics of a cosmic ray 
particle traversing the detector are the arrival direction 
(0,(j>), the particle identity and its kinetic energy or 
rigidity. These quantities were reconstructed combining 
the independent measurements provided by the various 
sub-detectors. The particle rigidity R = pc/Ze was 
given by the deflection of the particle trajectory 
reconstructed from hits in at least four out of the six 
tracker planes. The velocity j3 = v/c was determined 
from the transit time in the TOF along the track lenght. 
The particle charge \Z\ was obtained by the analysis of 
the multiple measurements of energy loss, as described 
in detail in the next section. 

In order to reject poorly reconstructed particles, a set 
of quality cuts was applied to the sample: 

1) particles were required to be measured by fully ef- 
ficient TOF scintillator paddles; this rejects events 
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passing through a single TOF paddle providing 
unreliable timing information 

2) tracks fitted with high \ 2 were removed according 
to a rigidity dependent cut-off value on the x 2 ; 

3) tracks with poor half rigidity agreement were also 
removed, according to rigidity-dependent criteria; 

4) consistency between independent (3 and R mea- 
surements was required, being f3 (particle velocity) 
and R (rigidity) obtained by two distinct subdetec- 
tors; 

5) the absence of hits was required close to the track 
extrapolation on missing layers, namely on tracker 
layers with no hit belonging to the reconstructed 
track. 

For the present work we considered data collected during 
these two periods: (a) MIR-docking phase, 4 days while 
the shuttle was docked with the MIR space station, 
and (b) post-docking phase, 2 days pointing at fixed 
directions (0°, 20° and 45° with respect to the zenith). 

Data collected while AMS-01 was passing in the 
region of the South Atlantic Anomaly were excluded 
(latitude: -47° -=--4°, longitude: -85° -j-— 5°) and only 
downward-going particles (in the AMS-01 reference 
frame) were considered, within a restricted acceptance 
of 32°. Approximately 75,000 nuclei from Lithium to 
Oxygen were selected with the above mentioned criteria. 

IV. Particle identification 

A nuclear species is defined by its charge Z. In this 
analysis the charge of each species is determined by the 
energy losses recorded in the silicon layers, due to the 
limited dynamical response of the TOF scintillators. The 
ionization energy generated by a charged particle in a 
silicon microstrip detector is collected by a cluster of 
adjacent strips. Tracker clusters were recognized online 
and then re-processed with the reconstruction software. 
A multi-step procedure of normalization of the cluster 
amplitudes was performed. The method accounts for 
saturation effects, electronics response, particle inclina- 
tion and velocity dependence of the energy loss. The 
charge identification algorithm, applied to the corrected 
signals, was based on the maximum likelihood method 
which determines the best Z value corresponding to the 
maximum value of the log-likelihood function: 

k=6 

L(Z) = log 10 {l[Pl(x k ,l3)} (1) 

k=l 

The k-index runs over the silicon layers, x^ are the 
corrected cluster amplitudes as observed on k-th layer. 
The functions have to be viewed as the 

probabilities of a given charge Z with velocity (3 to 
produce a signal xt on the k-th layers. These probability 
density functions were estimated from a clean reference 

'For a given track, the half rigidities i?i and R2 are two different 
measurements of the particle rigidity; R± (R2) was obtained with hits 
from the upper (lower) silicon planes. 



sample of flight data. For a fixed j3, they have gaussian- 
like shape; their dependence on f3 is shown on Figf2] 
Due to inefficiencies in charge collection, some energy 
losses may produce charge responses which do not 
carry reliable information on the particle charge. For 
this reason, not all the six layer clusters were used to 
determine Z. Firstly, a selection of the tracker clusters 
was done according to (1) reasonable strip occupancy 
levels (clusters containing noisy or dead channels were 
removed) and (2) cluster morphology (double-peaked 
clusters and single-strip clusters were removed). Finally 
the best set of three clusters given by ft = {k\, ki, k^} 
was recognized, was provided by the maximum like- 
lihood method as a parameter as well as Z. Hence, the 
actual log-likelihood function used was: 

L(Z,fi) = log 10 {P| 1 -P| 2 -P| 3 } (2) 

where Z ranges from 3 (lithium) to 8 (oxygen) and O 
runs over all the {fci, k-i, k^} 3-fold combinations of the 
tracker signals. It should be noted that, although the final 
Z estimation was provided by three tracker clusters, all 
the selected clusters were processed in Eq|2] 
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Fig. 2: Signal amplitude (ADC) of the mean energy 
loss in the silicon tracker vs velocity. Nuclear families 
fall into distinct charge bands. MPVs of the P^(xk,f3) 
functions are superimposed for Z=3 to 8 (dashed lines). 

V. Charge contamination 

To obtain the differential fluxes impinging the de- 
tector, the influence of charge contamination must be 
accounted for. Each nucleus of charge 3 < Z < 8 
produces a charge estimation Z which is related to its 
true impinging charge Z by using a set of coefficients 
F*. 

The 6x6 matrix \\F\\ is diagonally-dominated, and 
each off-diagonal element represents the probability 
of a nuclear species Z to be mis-identified as Z due to 
interactions in the detector material and fluctuations of 
the energy loss: 

F§ = P(Z\Z) (3) 

Two different contributions produce a charge 
migration Z — > Z : 
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Fig. 3: Charge histograms in the B and C region. The 
signal amplitudes of Fig f2] are here equalized to /3 = 1. 
Boron and Carbon families fall in distinct charge peaks. 

1) after interacting in TOF material, an incoming 
nucleus Z may fragment and physically turn into 
Z. These events were typically removed by the 
anticoincidence veto. However, a fraction of them 
produces a clean track on the tracker, passing trig- 
ger and selection. Since Z<Z, the corresponding 
matrix is triangular. 

2) due to energy loss fluctuations and broadening 
of the charge response, the overlapping of the 
measured dE/dx distributions produces an 
additional Z Z migration probability. This 
effect is tipically symmetric and in most cases 
Z = Z±l. 

Both the effects were studied with a Monte Carlo sim- 
ulation where nuclear interactions and detector responses 
were reproduced. An upper limit to the contribution (2) 
was also estimated by the data. In principle, by the exact 
knowledge of F z coefficients, it would be possible to 
correct the measured abundances Z to get the corrected 
ones (e.g. inverting the ||.F|| matrix). Since the overall 
distortion is estimated to be 2% -j- 6% for the ratio B/C 
(i.e. smaller than the statistical uncertainties), this effect 
was included in the systematic errors. 

VI. Flux calculation 

Detector counts are relate to the differential energy 
flux as a function of the kinetic energy per nucleon E kn 
according to the following relation: 

ftZfr \ N z (E kn ) 

9 ^ kn} A z (E kn )-T z (E kn )-AE kn W 

where N z (E kn ) is the number of the particles detected 
after selection in the energy bin around E kn , A z (E kn ) 
is the detector acceptance in the selected angular range, 



T z (E kn ) is the exposure time and AE kn is the width 
of the energy window around E kn . 

The acceptance A z (E kn ) is the convolution of the ge- 
ometrical factor with the energy dependent efficiency e, 
which includes (1) trigger efficiency, (2) reconstruction 
efficiency and (3) selection efficiency. A Monte Carlo 
estimation of the acceptance was performed according 
to the formula Q: 

A {Ekn) ~ A n§^T) ( ) 

where NQ(E kn ) is the number of Z-charged nuclei 
generated in the energy bin E kn and N^(E kn ) the 
number of those selected. Trajectories were generated 
from an ideal plane placed above the detector, and A = 
47.78 m 2 sr is the top plane geometrical acceptance. 
To calculate A B and A c , 2 • 10 8 boron and carbon 
trajectories were generated in 2-^250 GeV/c momentum 
range, according to an isotropic distribution and a power 
law momentum distribution with spectral index 7 = — 1. 

The geomagnetic cut-off, varying in the 1-15 GV 
interval along the AMS-01 orbit, introduces different 
distortions of the measured energy spectrum for different 
nuclear species. The vertical cut-off was calculated in the 
eccentric dipole approximation (4|: 

R vc = 15.0cos 4 (A)/i? 2 (6) 

where A is the geomagnetic latitude and R is the 
distance from the dipole center in terrestrial radii units. 
To remove any distortion due to geomagnetic field, 
only particles with rigidity greater than 1.2R VC were 
accepted. This means that, for a given species Z in the 
energy window AE kn , only some geomagnetic regions 
were accessible by AMS-01 during the orbit. The 
exposure time T z (E kn ) is then defined as the overall 
time spent by AMS-01 in the regions above cut-off for 
a certain species Z in the energy window around E kn . 
These T z were also corrected for the trigger deadtime 
fraction. 

The B/C ratio is then derived as the ratio of the 
differential fluxes 

B/C = 1> B (E kn )/<l> c {E kn ) (7) 

Two distinct measurements were made: (a) in the 0.4 -f- 
2.0 GeV/n range, the E kn was directly obtained from f3 
measurement; (b) between 2.0 GeV/n and 19.0 GeV/n 
the E kn was calculated by the rigidity from tracker. 
In the latter case, since the conversion R — > E kn 
requires the particle mass, an assumption on the B 
and C isotopic compositions is needecQ According to 
various measurements and predictions (e.g. [5|, |6|, [8 |), 
the cosmic ray boron is composed of a non-negligible 

2 Rigidity is the quantity directly measured by the tracker; it was 
translated into kinetic energy per nucleon according to the relation: 
Ekn(R) = yjR 2 {Z/A) 2 + A/2 - M p , being M v = nucleon mass 
and A= mass number. 
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fraction of 10 B. We considered B as a 10 B+ 11 B isotopic 
mixture and assumed a 10 B isotopic fraction as 



Y B = 10 B/( 10 B + 11 B) = 0.35 ± 0.15 



(8) 



whereas the carbon flux is assumed to be composed 
of pure 12 C. The uncertainty on Yg is consequently 
transposed as additional errors on the B/C ratio (3%-j-5% 
relative errors). 

VII. Results and discussion 

The results for B/C ratio are presented in FigH] The 
ratio B/C is calculated in ten energy intervals between 
0.4 and 19.0 GeV/n. The energy range was decided 
according to restrictions due to signal saturations (at 
lower energies) and requirement of reasonable statistics 
(at higher energies). The binning is uniform in logarithm 
scale; for each data point the reference value (Ekn) is 
defined as the geometric mean of the bin margins. The 
results are compared with HEAO-3-C2 data Q collected 
from October 1979 and June 1980. Both the detectors 
operated in space during periods of similar solar activity 
and polarity; altitude and inclination of the orbit were 
also comparable. The two measurements are in sub- 
stantial agreement. A theoretical prediction, calculated 
with the GALPROP numerical software (conventional 
diffusion model modulated with $ = 450MV fW&) is 
also reported. The dashed line is referred to a diffusion 
model (D. Maurin et al. Q) where the acceleration and 
diffusion processes are considered in terms of a rigidity- 
dependent diffusion coefficient 

K ^= K »e{-£v) (9) 

and a power-law distribution R~ a is assumed at the 
source; 6 = 0.6 and a = 2.2 are used; the solar 
modulation is described in the force-field approximation 
using $ = 500MV, consistent with the period of the 
AMS-01 flight. 

Both statistical and systematic uncertainties are repre- 
sented in the Fig. 4. Most of the systematic uncertainties 
arising from the detector response cancel into the B-C 
ratio. Differences in the trigger and selection efficiencies 
of the two species are expected from the Z dependence 
of the delta-rays production and fragmentation effects in 
the detector material. Deep investigation on these effects 
has been undertaken with GEANT4 and FLUKA fTT] 
based simulations with the aim of obtaining additional 
estimations of acceptances and their related systematics. 

Finally, since a considerable amount of light nuclei 
were recognized in the same energy range, one may 
expect to extend this analysis to other nuclear species. 
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Fig. 4: B/C ratio from this work (filled circles), from 
HEAO-3-C2 experiment (open diamonds) and from two 
diffusion models (solid and dashed lines). 
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